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Abstract 

The cross section of hard diffractive dissociation of the pion into two jets 
is calculated. It is obtained that the distribution of longitudinal momenta for 
jets is not simply proportional to the profile of the pion wave function, but 
depends on it in a complicated way. In particular, it is shown that, under 
the conditions of the E791 experiment, the momentum distribution of jets is 
similar in its shape for the asymptotic and CZ wave functions, and even the 
ratio of the differential cross sections is not far from unity. 

We argue therefore that, unfortunately, the E791 experiment has not yet 
measured the profile of the pion wave function. For this, the experimental 
accuracy has to be increased essentially. 



1. The E791 experiment at Fermilab has recently measured 

the cross section of the hard diffractive dissociation of the pion into two jets. 
In particular, the distribution of the total pion longitudinal momentum into 
fractions y± and y 2 , (yi + y 2 ) = 1, between jets has been measured. The 
main purpose was to obtain in this way the information about the leading 
twist pion wave function (j) n (xi, x 2 ), which describes the distribition of quarks 
inside the pion in the longitudinal momentum fractions X\ and x 2 , {X1+X2) = 
1. 

The hope was based on the theoretical calculations of this cross section in 
[2]- [4]. It has been obtained in all these papers that the cross section is simply 
proportional to the pion wave function squared: da/dyi ~ {^(.Vi)^ ■ In such 
it would be sufficient to measure only the gross features of da/dy 
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to reveal the main characteristic properties of 4> n (x), and to discriminate 
between various available models of <p n (x). 

The purpose of this paper is to show that this is not the case. The real 
situation is much more complicated, with da/dy depending on ^(x) in a 
highly nontrivial way. We give below (in a short form) the results of our 
calculation of this cross section. 

2. The kinematics of the process is shown in fig. [TJ 

We take the nucleon as a target, and the initial and final nucleons are 
substituted by two soft gluons with momenta q\ and q 2 . The lower blob in 
fig. 1 represents the skewed gluon distribution of the nucleon, g^{u). 

The final quarks are on shell, carry the fractions y\ and y 2 of the ini- 
tial pion momentum, and their transverse momenta are: (kj_ + (qj_/2)) and 
(— + (q_i_/2)), q± <C k±, where q± is the small final transverse momentum 
of the target, while k± is large. 

The upper blob M in fig.l represents the hard kernel of the amplitude 
which includes all hard propagators. For calculation of M in the leading 
twist approximation and in the lowest order in a s , the massless pion can be 
substituted in all diagrams by two massless on shell quarks with the collinear 
momenta Xip n and x 2 V-n and with zero transverse momenta, as account of 
primordial virtualities and transverse momenta results only in higher twist 
corrections to M. The leading twist pion wave function cj) n (x,fj, ) describes 
the distribution of these quarks in momentum fractions X\ and x 2 . [] 

The hard kernel M is proportional to the scattering amplitude of two 
initial collinear and on shell quarks of the pion on the on shell gluon: 

d(xxp n ) + u(x 2 p n ) + g(qi) -> dipt) + u(p 2 ) + gigz) 

In lowest order in a s (k±) M consists of 31 connected Born diagrams, each 
one is ~ 0(a^(k±)) and contains exactly three hard propagators (except for 
one diagram with the 4-gluon vertex, which has only two). Two diagrams 
are represented explicitly in fig. and fig. || 

In the c.m.s. and to the leading twist accuracy, the initial and final 
soft gluons can be considered to be on shell, with transverse polarizations, 
carrying fractions (u + £) and (u — £) of the mean nucleon momentum P. 

1 As usual, on account of leading logs from loops the soft pion wave function (f) w (x, /x D ) 
evolves to ff (x, n) ~ J M d 2 Z^4 , T (x, l±), where fx is the characteristic scale of the process. 
And the same for the gluon distribution: g^{u, fJ, a ) — > g^(u,/j,). 
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Figure 1: Kinematics and notations 



f\ The nucleon can be considered as being massless, spinless, and its skewed 
gluon distribution g^(u,t,fi ) = g^(—u,t,fi ), —1<u<1, is defined as f\ : 



(pK'>i)^> 2 )i^V 



5 ab r 1 du 



-9xu- 



■ x 



x I ^ e -i(u+0(Pvi)+i{u-t)(Pv 2 ) _|_ e -i(«+0(-P«2)-N(«-0(-P«l)^ . 

This is equivalent to the standard definition (see [5-7]): 



(1) 



1 (e -* 



+€)(-Pfl)+*(«-C)(^2) I e -i(u+S)(- P «2)+i(M-S)(PDl) 



The kinematical variables are defined as: 

?i = (« + 0?, ?2 = («-0^ P = (P + P')/2 



(2) 



2 The small skewedness, £ <C 1, is always implied. It is typically: £ ~ 10 -2 , in the 
Fcrmilab experiment. 

3 Here and below the gluon string is always implied in gauge invariant definitions of 
bilocal operators, both for the gluon distribution of the nucleon g^{u), and for the pion 
wave function (j>„(x). 
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A = [qi - q 2 ) = 2 f P, £ = ^ , *l = , ^2 



2yms 



2£ 



2£ 



(3) 



*i - z 2 = 1, 2 (p ff A) = M 2 



2/12/2 



According to the well developed approach to description of hard exclusive 
processes in QCD [8]- [11] (see [12] for a review), all hard gluon and quark 
lines in all diagrams (see e.g. figs. 2 and 3) have to be written down explicitly 
and substituted by their perturbative propagators. In other words, the hard 
momentum flow have to be made completely explicit and these hard lines 
of diagrams constitute the hard kernel M. They should not be hidden as (a 
derivatives of) "the tails" of the unintegrated pion wave function ^^(x, l±), 
or of the "unintegrated gluon distribution". This is, first of all, what differs 
our approach from previous calculations of this process in [2-4] where, besides, 
most of diagrams were either ignored or calculated erroneously. 

Due to this the function *$? w (x, l±) ( and the unintegrated gluon distribu- 
tion) never appears explicitly in our calculations. It always enters the answer 
for each of 31 Born diagrams implicitly only through the function <p n (x,fi), 
i.e. only in the standard integrated form: (j) n (x, //) ~ J M d 2 l^ l±). Q 

So, the structure of the amplitude is (symbolically): 

T ~ (P'lA 1 ■ A ± \P) ® (^M^ 2 ) ® (0|w • d\n-), (4) 

where the first matrix element introduces the skewed gluon distribution of 
the nucleon g^(u), ipi and tp 2 are the free spinors of final quarks, "M" is the 

4 While a small number out the whole set of 31 diagrams can possibly be reinterpreted 
through "the tails" of the unintegrated pion wave function ^ v (x,l±) or of the "uninte- 
grated gluon distribution" , this is definitely not the case for most of diagrams. So, we see 
no much meaning to proceed in this way. 

In a sense, the contribution of each diagram is essential in obtaining the numerical 
description for the "y" -distribution of jets in sect. 4. In our approach all 31 diagrams are 
treated on equal footing. 

5 In H the authors tried to use the evolution equation for the pion wave function to 
obtain its "tail". For instance, for the asymptotic wave function it was obtained in this 
way: ^^(x, k±) ~ (^^(x, [i)/k±. It remains unclear for us how it is possible to obtain 
such a result from the evolution equation for the asymptotic wave function which looks 
as: d4Z?Hx, f i)/dln f i = J dyV(x,y)r/ v (y^) = 0- 
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hard kernel, i.e. the product of all vertices and hard propagators, the last 
matrix element introduces the pion wave function 4> n (x), and (g> means the 
appropriate convolution. 

As an example, let us consider the diagram in fig. 0. Proceeding in the 
above described way (see the appendix), one obtains the contribution to the 
amplitude (the Feynman gauge is used for the hard gluon): 

T = -— — f 1 dXl f 1 — 9 ^ U>> (b) 

9 y 2 Jo x x x 2 J-i (u — 0i u + ' 



^ = ^^^A(^A75^)^^, A = A, 7 „ (6) 

where ip 1 and ip 2 are the free spinors of the final quarks, E sp j ns |^ 1 A7 5 ip 2 \ 2 = 
2k 4 L /(yiy 2 ), "ij" are their colour indices, and /„. ~ 130 Me V is the pion 
decay constant. 

As it is expected that, for this process, the imaginary part of the am- 
plitude is the main one at high energy, we show explicitly in ch.3 only its 
value. For the diagram in fig. [| this gives (in all diagrams the terms ie are 
introduced into denominators through s — > s + ie, i.e. £ — > £ — ie): 

I6 2irsuj y 1 f 1 dx x (p n (x) 

As a final example, let us also consider the diagram in fig. as it gives 
(together with the mirror diagram obtained by q\ «-> q 2 ) the enhanced con- 
tribution [] : 

T ^ = _Uo_ f 1 dxi fajx) r 1 dug^ju) N 

3 ViV2JQ xxx 2 J-i(u-£)(u + £,)[zi(x 1 -yi)-x 1 y 2 \ 



N = [-8 ziz 2 ] + (8 2/i2/2 - 3 - x x y x - x 2 y 2 ) + (zi + z 2 )(x l - y x ) , (9) 



ImT 3 = 8- r / : — — 6(|xi — 2/1 > o) -\ , (10) 

kj_ Jo x\x 2 \x\ — 2/1 1 



6 The analytic form of the enhancement depends on the behaviour of the gluon distribu- 
tion g^(u), see eq.(10). With g^(u) ~ const the enhancement is logariphmic, ~ ln(s/fcj_), 
but with g^(u) from eq.(18) it is only numerical. 
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u = d XlV2+Xm ), S = kl/s, 



(ii) 



where only the enhanced term is shown explicitly in eq.(10), which originates 
from the term in square brackets in eq.(9). [] 

3. Proceeding in the above described way and summing up the 

contributions of all 31 Born diagrams, one obtains for the cross section: 



do 



N 



T\ 



dy 



8 (2tt)5 S 2 i i ym 



d k±d q±, T~iImT = i- 



2ir sco. 



£ ft(0«. (12) 



n 



dxi ^{xi) (Si 



E 3 + S 



, 



(13) 



Si 



9iW 



xix 2 \x 1 - y ± \ gt(£) 



0(|xi - yi\ > 6) 



+ (yi 2/2), 



(14) 



ac? z§ l/i I/2 



-(xix 2 + ym)+ 



+ 



Lr 



1 - yi|(xi - y 2 ) 2 9j ^r 9(|xi - yi| > 5) + (3/1 <-> y 2 ) 



(15) 



9 x\x\y x y 2 



Fi - 2/1 Pttt ©( Fi - 2/i I > 0) + (2/1 2/2) 

&(0 



,(16) 
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S 4 = — 



1 



9 xix 2 ym 



dg^{u) I du\ u= 

#(0 



(17) 



The expressions (12)-(17) constitute the main result of this paper. 



7 As the 8 - function in eq.(10) excludes the region of x\ too close to y±, the integral 
is convergent and only some enhancement remains (see the footnote 6). 

At the conditions of the E791 experiment (and even at much larger energies, see fig. 6), 
this enhancement is by no means sufficient to neglect contributions of all other diagrams. 
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Let us note that while the separate terms in / <ix0 7r (x)S 2 are logarith- 
mically divergent at Xi >2 — ^ 0, it is not difficult to see that the divergences 
cancel in the sum, so that the integral is finite. And the same is valid for 
S3. This is an important point as it shows that the whole approach is self- 
consistent, i.e. the hard kernel remains hard and the soft end point regions 
£1,2 — > give only power suppressed corrections. 

4. In this section we present some numerical estimates of the cross 

section, based on the above expressions (12)-(17). Our main purpose here 
is to trace the distribution of jets in longitudinal momentum fractions yi,y2 
depending on the profile of the pion wave function (j) n (x). 

a) As the calculations were performed in the leading twist approximation 
which becomes applicable at sufficiently large k± only, we take kj_ = 2 GeV. 
This assumes that the higher twist effects are not of great importance at such 
a value of k±, and gives a possibility of comparison with the E791-data. 

b) For the skewed gluon distribution g^(u,t,fi) of the nucleon at t — 
—qj_ ~ we use the simple form (as we need it at \u\ > £ only, and because 
9t;( u ) 9o(u) at \u\ > £): 

gt(u, t = 0, n ~ k± ~ 2 GeV) | u > e ~ u -°- 3 (l - uf. (18) 

This form agrees numerically reasonably well with the ordinary, g (u, /i ~ 
2 GeV), and skewed, g%(u,t = 0,/i ~ 2 GeV), gluon distributions of the 
nucleon calculated in [13] and [14] respectively (in the typical region of the 
E791 experiment: \u\ > £ ~ 10 -2 ). 

The detailed consideration of nuclear effects is out the scope of this paper. 
So, we simply assume that these effects result mainly in an overall factor (see 
1): 

da A (t ~ -q\) ~ do N (t = 0) \AF A {t)\ 2 , (19) 

where F A (t) = exp{bt/2}, b = (R A )/3, is the nuclear form factor. 

c) As for the pion leading twist wave function, (p n (x, fi), we compare two 
model forms: the asymptotic form, (p^f y (x, fi) = 6x1X2, and the CZ-model 

15| . The latter has the form: (p^ z (x, \i Q ~ 0.5 GeV") = 30xiX2(xi — X2) 2 , 
at the low normalization point. Being evolved up to the characteristic scale 
ji ~ k±_ ~ 2 GeV, it looks as (see fig. f|): 

(f/{x, fi-2 GeV) = 15 xix 2 f(xi - x 2 f + 0.2l . (20) 



As the final quarks are free, the wave function of this two quark system is the asymp- 
totic one, and it does not evolve in the leading log approximation. 
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Figure 4: Profiles of the pion wave functions: a) (f>^ (x,/jl ~ 0.5 GeV) = 
30xix 2 (xi - x 2 ) 2 - dotted line; b) <p^ z (x,fi ~ 2Gey) = 15xix 2 [0.2 + (xi - x 2 ) 2 ] 
- solid line; c) </>" SJ/ (x) = 6x1X2 - dashed line. 
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The results of these numerical calculations are then compared with the 
E791-data, see fig. |[ 

It is seen that, unfortunately, while two pion wave functions are quite dif- 
ferent, the resulting distributions of jets in longitudinal momenta are similar 
and, it seems, the present experimental accuracy is insufficient to distinguish 
clearly between them. Moreover, even the ratio of the differential cross sec- 
tions is not much different from unity: da asy /da cz ~ 1.2 at y\ = 0.5, and 
the same ratio is ~ 0.7 at y\ = 0.25. 

In such an unhappy situation, the theoretical calculations should be also 
performed with a maximal possible accuracy (to account additionally for: the 
quark distributions, higher twist corrections, hard loop corrections, nuclear 
effects etc.). 

We also show in fig. |6| the same distributions with the pion energy ten 
times larger, E n = 5TeV, k± = 2 GeV. It is seen that even this does not help 
much (as the form of the distribution weakly depends on the pion energy). 
The same ratios of the cross sections are here ~ 1.7 and ~ 0.7 respectively. 

Recently the Coulomb contribution to the cross section has been calcu- 
lated in p|. PI Its value for Pt is: da e ^ ctr / dk 2 ± dy < 10~ 5 mbarn ■ GeV~ 2 , 
for E n = 500 GeV, k± = 2 GeV. Using the above given formulae, one ob- 
tains for the strong cross section at the middle point y = 0.5: dapf jdk\dy ~ 
2.5- 10~ 2 mbarn- GeV~ 2 with the same parameters, and ~ 0.25 mbarn- GeV~ 2 
at E n = 5TeV. It is seen that the electromagnetic contribution is small. 



Note added: After this work has been completed, the paper |T7J on 
the same subject appeared. Comparison shows that, although the qualitative 
conclusions are similar, the analytic expressions for the scattering amplitude 
differ essentially in this paper and in fl7j . 
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Figure 5: The y-distribution of jets calculated for k± = 2GeV, = 500 GeV and 
with the pion wave functions: 4>^ z (x, n ~ 2GeV) - solid line, 4>® sy (x) - dashed 
line. The overall normalization is arbitrary, but the relative normalization of two 
curves is as calculated. The data points are from the E791 experiment [1]. 
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Figure 6: The same as in Fig.5, but with E w = 5TeV. 
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Appendix 



The purpose of this appendix is to give (somewhat schematically) some 
details of how the diagrams have been calculated. As an example, let us 
consider the simplest diagram in fig. 2. The momentum of the virtual gluon 
is: k = (xip n — pi). The momenta of virtual quarks are: cr = (— X2p w — Z]A) 
and p = (—X2p n — A), with g 12 = £i )2 A, see eq.(3). The overall denominator 
is therefore: D = k 2 a 2 p 2 = (—XiX 2 .y2Zi)(2p 7T A) 3 . This diagram looks (in the 
operator form) as: 



To 



. (47ra s . 
D 



—X a 



At A: 



(21) 



The matrix element of this operator between the initial and final states is 
factorized then into three parts. The first part is the matrix element of the 
gluon fields between the initial and final nucleons, this introduces the gluon 
distribution, see eq.(l): 



\bc 



(22) 



(P'\Ai' b A^ C \P) - . aw g A{u _ 0{u + - 
The second part is the matrix element of two quark fields d a uL between 



the pion and vacuum, this introduces the pion wave function Q : 



(014410 



)a/3 



(23) 



p\ i 3 4 

The last part is the matrix element of two remaining quark fields between 
the vacuum and the final state of two free quarks. This introduces the Dirac 
spinors ip 1 and f/V The result is integrated then over "x" and "u". So, one 
obtains: 



dx\ I du (47ra; s ) 



9d 



u 



12 32(u-0(w + V3 



(^iM^ 2 



1 

D 



(2^A) 
D 



8ij (^M^ 2 ) ,(24) 
(^ 1 A 75 ^ 2 ) • (25) 



On the whole, one obtains the eqs.(5,6). All other 30 diagrams have been 
calculated in the same way. 

11 Eq.(23) is a shorthand for a standard definition of the pion wave function <j> v (x) 
through the matrix element of the gauge invariant bilocal quark operator, see [8, 9], in the 
same way as eq.(22) is a shorthand for eq.(l). See also the footnote 3. 
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